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We show that, based on the idea of string fragmentation, the production rates of light flavored mesons and 
baryons originating from fragmentation can be described by the spin, the binding energy of the particle, and a 
strangeness suppression factor. Apart from a normalization factor, e+e — data at different center-of-mass energies 
can be described simultaneously. Applying to the heavy flavor production, we find that our predictions are in 
good agreement with data. 



1 Introduction 

The soft processes of hadronization are not cal- 
culable with a perturbative approach and instead 
rely currently on phcnomcnological models. The 
most successful of these models are the stringOl and 
the cluster B fragmentation models, implemented 
in the_Monte Carlo programs JETSET0 and HER- 
WIG respectively. However, these models ei- 
ther require a large number of free parameters in 
order to reproduce the measured hadron produc- 
tion rates ( JETSET) , or do not give a satisfactory 
description of baryon-data (HERWIG). A review 
may be found in Ref.Erfl. 

In this paper, we first show some regularities 
in the hadron production rates measured at LEP. 
Then based on the idea of string fragmentation, 
we deduce a simple formula to describe the LEP 
data. The formula is also applied to data obtained 
at center-of-mass energies around 10 GeV and 29- 
35 GeV, and to the heavy flavor production. 

2 Overview of LEP data 

Thanks to excellent performance of the detectors 
and high statistics available, very careful work by 
all four LEP experiments has given a very com- 
plete picture of the production of identified parti- 
cles from e + e _ annihilation. The measured pro- 
duction rates per hadmnic Z event for the identi- 
fied particles at LEP Era are listed in Table |[ 

Studies of general features of particle produc- 
tion, such as the fraction V/(V + P) of mesons 
produced in spin-1 states and the strangeness sup- 
pression factor 7 S = s/u, provide useful informa- 
tion about the fragmentation process. The ratio 
V/ (V + P) for primary mesons is expected to be 
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equal to 0.75 from simple spin counting. From 
the rates and the primary fractions given in Ta- 
ble g and Ref. 0, we obtain a value of 0.46 ± 0.04, 
0.42 ± 0.03, 0.56 ± 0.04 and 0.75 ± 0.04 for u(d)-, 
s-, c- and b- mesons respectively. The low value 
of V/(V + P) for light flavored mesons could be 
explained by mass differences between the vector 
and pseudoscalar mesons, i.e. by the relatively 
larger binding energy of pseudoscalar mesons. For 
c-mesons the measured ratio of V/ (V + P) is also 
low, while for b-mesons it agrees well with the ex- 
pected value of 0.75. We will discuss this later in 
more detail. 

Figure 1 shows the production rates of pri- 
mary hadrons (measured rate x fraction from 
JETSET as listed in Table 0), divided by the 
spin factor (2J+1), as a function of the hadron 
mass. We see similar behavior ^] (steps) for 
hadrons belonging to the same multiplet, such 
as (p/u,K*,<p), (p,A,E,S) J (A,£*,S*,Q) and 
(f2,K2,f^). The mass difference of hadron pairs 
which are in the same multiplet but differ by one 
in the strangeness is in the range of 0.1-0.15 GeV 
(except m K —m n ). This is close to the mass differ- 
ence between the s and the u(d) quark, showing 
that the binding energy of hadrons in each pair is 
about the same. Unlike the case of V/(V + P), one 
could expect the 7 S value determined from differ- 
ent hadron pairs is about the same. 

In the string fragmentation model, one ex- 
pects the strangeness suppression factor 7 S to be 
around 0.3. This parameter can be measured 
from the production rates of strange compared 
with non-strange hadrons, and from the momen- 



a From Fig. 1 one can see that it is not possible to de- 
scribe all the data with just an exponential function of the 
hadron mass (or mass squared). 
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Table 1: Average particle production rates (excluding 
charge conjugates and antiparticles), compared to the cal- 
culated values (see next section). The fraction of primary 
hadrons obtained from the fit and JETSET is also shown. 



Particle 


Rate 


Rate 


Prim. 


Prim. 




Measured 


Calc. 


Frac. 


Frac. 








Calc. 


JETSET 


7T U 


9.19 ± 0.73 


9.77 


0.16 


0.14 


7T + 


8.53 ± 0.22 


8.70 


0.18 


0.16 


K° 


1.006 ± 0.017 


1.008 


0.25 


0.30 


K + 


1.185 ± 0.065 


1.100 


0.23 


0.27 


n 


0.95 ± 0.11 


0.85 


0.30 


0.33 


n' 


0.22 ± 0.07 


0.11 


0.59 


0.79 


to (980) 


0.140 ± 0.034 


0.080 


0.99 


0.93 


— n 

P 


1.29 ± 0.13 


1.12 


0.47 


0.54 


K 


0.380 1 0.021 


0.390 


0.49 


0.60 


K ^ 


0.358 ± 0.034 


0.395 


0.50 


0.60 




1.11 ± 0.14 


1.05 


0.48 


0.57 


<p 


0.107 ± 0.009 


0.092 


0.64 


0.70 


f 2 (1270) 


0.25 ± 0.08 


0.19 


0.78 


0.96 


KJ(1430)° 


0.095 ± 0.035 


0.051 


1.00 


0.98 


£,(1525) 


0.020 ± 0.008 


0.018 


1.00 


0.98 


P 


0.49 ± 0.05 


0.54 


0.12 


0.56 


A 


0.186 ± 0.008 


0.165 


0.12 


0.44 


E° 


0.0355 ± 0.0065 


0.0387 


0.40 


0.86 


E+ 


0.044 ± 0.006 


0.037 


0.42 


0.88 




0.0129 ± 0.0007 


0.0118 


0.41 


0.75 


A++ 


0.064 ± 0.033 


0.069 


0.69 


0.95 


E(1385) + 


0.011 ± 0.002 


0.016 


0.91 


0.92 


H(1530)° 


0.0031 ± 0.0006 


0.0050 


0.94 


0.94 


Q _ 


0.0008 ± 0.0003 


0.0015 


0.88 


0.92 


A(1520) 


0.0107 ± 0.0014 


0.0129 


0.71 




turn spectrum of strange 


mesons 


The 


results 



are very consistent wit, 
may be found in Rcf. 



the expectation (a review 
). This suggests that the 
strangeness suppression occurs at the quark level. 
In the following we consider that the particle pro- 
duction proceeds in two stages, namely quark pair 
production in the color string field and successive 
recombination. 

3 Analysis 

Quark pair production in the color string field 
can be considered as a tunneling process. The 
probability to produce a qq pair is proportional 
to exp(— 7rm^/ft), where m q is the (constituent) 
quark mass, and k the string constant. We as- 
sume that the probability of quarks recombined 
to a hadron with the mass Mh is proportional to 
exp(—E bind /T), where E bin d = Mh~Y,i m <a is tne 
hadron binding energy |, and T the effective tem- 
perature in hadronization. The production rates 
of light flavored mesons and baryons from frag- 
mentation can be described as 



<iv>=c.H^±i.( 7s) ^ 
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Figure 1: Production rates of primary light flavored 
hadrons at LEP energy, divided by the spin factor (2J+1), 
as a function of the hadron mass. 

where 7 s =exp(— Tr{m 2 s — m„)//c) is the strangeness 
suppression factor, N s the number of strange 
quarks contained in the hadron, and J the spin 
of the hadron. C is an overall normalization fac- 
tor, and Cb is the relative normalization factor 
between mesons and baryons (for mesons Cs = l). 
Equation |l| can also be applied to mixed states of 
the SU(3) octet and singlet, such as 77 and 7/, by 
adding up the uu(dd) and ss contributions. For 
this purpose we use the mixing formulae and an- 
gles given in Ref. M . 

Total hadron rates are calculated as follows. 
At first the number of light flavored hadrons pro- 
duced from fragmentation is calculated by using 
Eq.[l]. For hadrons which contain a primary quark 
q, we use Eq.[l] to determine their relative ratios, 
and then get their rates by normalizing the sum of 
the rates to the qq fraction, Tqq/Thad, which can 
be calculated by the Standard Model. All light 
flavored hadrons up to a mass of 2.5 GeV ia the 
meson and baryon summary table of Ref. liil arc 
included in the calculation. In the next step we 
let all these primary hadrons decay according to 
thcir_decay channels and branching ratios given in 
Rcf.Eil. The decay chain stops when fj,, n, K ± , K° L 
or stable particles are reached. 

In the fit we choose j s . Am = m s — m u , T, C 
and Cb as free parameters^ (we assume m u —md). 



b It can be ascribed to the hyperfine interaction 0. 



c Since cxp(-E bind /T) = exp(^\ m u /T) ■ exp(-(A/ h - 

^2ii^ m ii ~ mn ))/T)> the factor exp(y"\ m u /T) can be ab- 
sorbed in C and Cb ■ The error function of the fit is mainly 
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Table 2: Results of the fit to LEP data and to data obtained 



at different center-of-mass energies. 



Parameters 


sfi = 91 GeV 


Simultaneous Fit 


7s 


0.29 ± 0.03 


0.29 ± 0.02 


Am (GeV) 


0.150 ± 0.029 


0.161 ± 0.024 


T (GeV) 


0.289 ± 0.020 


0.298 ± 0.015 


C 


0.209 ± 0.041 


Cgi = 0.218 ± 0.034 * 


C B 


10.0 ± 1.0 


11.0 ± 0.9 


x'/dof 


56.2/19 


155.8/57 



* C 30 = 0.124 ± 0.020, Cio = 0.049 ± 0.008 



The results of the fit to the LEP data]] for a typical 
s quark mass of 0.5 GeV are listed in Table ^. The 
calculated total rate and primary fraction for dif- 
ferent hadrons are listed in Table [|. We see a good 
agreement between the measured and calculated 
rates except fbr_some decuplet baryons. As men- 
tioned in Ref.lZrtil, experimental errors for the decu- 
plet baryons are still large and there are discrep- 
ancies between experiments. The large % 2 value 
of the fit is mainly due to the decuplet baryons. 
If they are excluded in the fit, the x 2 /dof is then 
equal to 22.6/15 while the fit results remain essen- 
tially unchanged. However, the difference between 
data and our predictions for the decuplet baryons 
might suggest that baryon production is not de- 
scribed properly in our approach, since baryons 
are considerably more complicated objects than 
mesons. If diquark production!!! is the main mech- 
anism for baryon production, it is then more suit- 
able to use the effective diquark masses than the 
quark masses in Eq|j]. This will, however, increase 
considerably the number of free parameters. 

We also perform a simultaneous fit to data 
at different center-of-mass energies with 7 free pa- 
rameters: 7 S , Am, T, Cb, Cm, C30 and Cio, where 
C91, C30 and Cio are the overall normalization 
factor at the center-of-mass energy 91, 29-35 and 
10 GeV, respectively. The fit results are listed in 
Table ^ and shown in Fig. 2. The value of j s , 
Am, T and Cb is very consistent with that ob- 
tained from the LEP data alone, showing that the 
value of these parameters is independent of the 
center-of-mass energy. 

The 7s value determined from the fit is in an 
excellent agreement with the expectation. The 
value of Am = m s — m u required from the 
fit is consistent with that obtained from hadron 
masses Ej and from QCD calculations 113. While 
the values of 7 S , Am and T do not depend on the 
choice of the m s value, the values of C and Cb are 

sensitive to the change in the mass difference m s — m u . 
d The new result for A(1520) is not included in the fit. 



• Vs=91 GeV 
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Figure 2: Production rates of light flavored hadrons orig- 
inating from fragmentation at different center-of-mass en- 
ergies, multiplied by the factor Cb /(is 1 " (2J + !))> as a 
function of the binding energy of hadrons. 

correlated with the m s value. However, the calcu- 
lated baryon rates do not depend on the choice of 
the m s value since they are determined by both m s 
and Cb- From our analysis we find that baryons 
contribute 10.9% of the primary produced hadrons 
at LEP energy. 

Fractions of primary hadrons determined from 
the fit are listed in Table [l]. They are in general 
lower than those predicted by JETSET due to the 
inclusion of orbitally excited states up to a mass of 
2.5 GeV in our analysis. JETSET includes some 
of orbitally excited meson states (only those with 
L = 1), but no orbitally excited baryon states. 
From our analysis we obtain the fraction of pri- 
mary produced states with orbital excitation to 
be 27% for mesons and 39% for baryons at LEP 
energy. 

Using the T value obtained from the fit, the 
ratio V I (V + P) can be calculated according to the 
mass difference of vector and pseudoscalar mesons. 
We obtain a value of 0.38, 0.44, 0.65 and 0.72 for 
u(d)-, s-, c- and b-mesons respectively, which are 
in good agreement with the measurements. 

4 Predictions of heavy flavor composition 

Using the fit results obtained with light flavored 
hadrons, the relative production rate of heavy fla- 
vored hadrons can be predicted by our method. 
The absolute value of the rates can be deter- 
mined by normalizing the sum of the rates to the 
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Table 3: Predictions of average production rates for heavy 
flavored hadrons, compared with LEP data. 



Hadrons 


Prediction 


Measurement 


D u 


0.242 


0.221 ± 0.012 


D+ 


0.092 


0.087 ± 0.008 


D* + 


0.114 


0.088 ± 0.006 


D+ 


0.054(±0.005) a) 


0.041 ± 0.008 


A+ 


0.026 


0.037 ± 0.009 


B° 


0.091 


0.097 ± 0.026 


D*/(D* + D) b) 


0.66 


0.56 ± 0.4 


B*/(B* +B) 


0.70 


0.75 ± 0.04 


< D° + D%° > c Br(D 


*° -t D* + 7T-)/ < D*+ > c 6) 




0.124 


0.10 ± 0.03 


B:* d /B u , d 


0.39 


(0.24 ± 0.03)/ c) 


(Bi +B*)„, d /B„, d 


0.228 


0.216 ± 0.033 


(b s1 + B; 2 )/B;* d 


0.110 


0.142 ± 0.055 


(b s1 + b; 2 )/b+ 


0.086 


0.052 ± 0.016 


D"/c-jct 


0.593 


0.570 ± 0.046 


1 J l C J 1; L 


0.237 


0.249 ± 0.026 


rt*~f / f 

L) 1 C-JCt 


0.272 


0.241 ± 0.015 


L) s /c-jet 


0.101(±0.025) a) 


0.128 ± 0.027 


B°/b-jet 


0.108(±0.030) d) 


0.122 ± 0.031 


A+/c-jct 


0.069 


0.076 ± 0.044 


Ab/b-jct 


0.073 


0.076 ± 0.019 e) 


(£(, +S 6 *)/b-jot 


0.060 f) 


0.048 ± 0.016 


D s i/c-jet 


0.012 


0.016 ± 0.006 


(D sl +D; 2 )/c-jct 


0.028 




(B 8l + BJ 2 )/b-jet 
(Di + D') Utd /c-jet 


0.035 


0.037 ± 0.012 


0.170 


0.173 ± 0.053 


(Bi +Bj) Uid /b-jet 
D**/c-jot 


0.188 




0.38 




B**/b-jct 


0.38 




c-baryon/c-jet 


0.089 




b-baryon/b-jct 


0.091 


0.115 ± 0.040 



^ for Br(D" — D S X) = 0.5(±0.5) 
b ' excluding D. D* ^lid D** from B 

c) / = 1-2 (sec Ref. 0) 

d) for Br(B" -» B S X) = 0.5(±0.5) 

e) for Br(A 6 — A+TvX) = 0.1 
n - 0.052 if m s> ,. E » - m A , 



decays 



measured in Ref.Lfl arc used 



^cc/^had or T b i/T had , which can be calculated by 
the Standard Model. We use all heavy flavored 
hadrons in the JETSET 7.4 table, which includes 
orbitally excited mesons with L = 1, such as D** 
and B**, but no orbitally excited heavy flavored 
baryons (for more discussions see Ref.la). 

As an example we consider the heavy flavor 
composition at LEP energy. Our predictions _tia- 
gether with the corresponding measurements uO, 
are listed in Table |j. There is a good agreement 
between the predicted and measured values. The 
main uncertainty on the prediction is due to the 
limited knowledge on properties of orbitally ex- 
cited heavy flavored hadrons, such as mass, decay 
modes, and branching ratios Without taking 
into account the production of orbitally excited c- 
and b-baryons, the fractions c-baryon/c-jet and b- 
baryon/b-jet are calculated to be around 9%. The 
calculated values would be 30-40% higher □ if or- 
bitally excited c- and b-baryons could be included, 
which are however barely known. 



5 Conclusions 

We have shown that the production rates of light 
flavored mesons and baryons in e + e _ annihilation, 
which span a range of four orders of magnitude, 
can be described by a simple approach based on 
the idea of string fragmentation. Data at different 
center-of-mass energies can be described simulta- 
neously apart from a normalization factor, which 
reflects the rise of multiplicities with increasing en- 
ergy. From the fit we determine the strangeness 
suppression factor to be 0.29 ± 0.02. Applying 
to the heavy flavor production, we find that our 
predictions are in good agreement with data. A 
comparison of our approach with other recently 
proposed approaches can be found in Ref. Q. 

The proposed approach may provide insight 
into the hadron production mechanism. To further 
check the approach, precise data on the produc- 
tion of the decuplet baryons and orbitally excited 
hadrons are needed. In particular, with better un- 
derstanding of the baryon production mechanism, 
the description of baryon production in our ap- 
proach can be improved. 
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